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ABSTRACT 

We examine the Lya absorber population at z < 0.3 detected in spectra of the QSOs 
PG0953+415 and H1821+643 taken with the Space Telescope Imaging Spectrograph aboard the 
Hubble Space Telescope. We compare their statistical properties to those in carefully-constructed 
mock quasar spectra drawn from a cosmological hydrodynamic simulation of a ACDM universe. 
We find very good agreement in the column density and 6-parameter distributions, down to the 
smallest observable absorbers with TVhi ~ 10 123 cm~ 2 . The observed absorber population is com- 
plete for TVhi > 10 13 cm~ 2 , with a column density distribution slope of /3 = 2.04 ± 0.23 and a 
median 6-parameter of 21 km/s above this limit. The intergalactic gas giving rise to these weak 
absorbers is analogous to that at high redshift, located in diffuse large-scale structures that are 
highly photoionized by the metagalactic UV flux, though a greater number arise within shock- 
heated warm gas. The density, temperature, and column density of these absorbers follow similar 
relationships to those at high redshift, though with substantially larger scatter due to the shock- 
heated gas. The 6-parameters typically have a significant contribution from thermal broadening, 
which facilitates a measurement of the low-z IGM temperature as traced by Lya absorbers. From 
our simulation we estimate Tiqm ~ 5000 K, with an upper limit of 10 4 K, at the mean density. 
The agreement in predicted and observed amplitude of the column density distributions allows 
us to measure the H I photoionization rate at z = 0.17 to be Phi = 10 _13 - 3±0 - 7 s~ 1 (estimated 
modeling uncertainty), close to predictions based on quasar properties. 

Subject headings: cosmology: observations — cosmology: theory — intergalactic medium — quasars: 
absorption lines — quasars: individual (PG0953+415,H1821+643) 



1. Introduction 

Neutral hydrogen along the line of sight to 
distant quasars produces numerous Lyman al- 
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pha (Lya) absorption features in quasar spectra, 
known as the Lya "forest" (Lynds 1971). At 
redshifts z > 2.5, high-resolution optical spec- 
troscopy with HIRES (Vogt et al. 1994) on the 
Keck 10m telescope has enabled the statistics of 
Lya absorbers to be determined to high preci- 
sion (Hu et al. 1995; Lu et al. 1996; Kim et al. 
1997; Kirkman & Tytler 1997). Recent obser- 
vations using UVES at the VLT have probed ab- 
sorbers at redshifts 1.5 < z < 2.5 (Kim, Cristiani 
& D'Odorico 2000). In conjunction with these 
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observations, hydrodynamic simulations of the in- 
tergalactic medium (IGM) have forwarded a self- 
consistent physical picture for the origin of Lya 
forest absorbers (-/Vhi < 10 17 cm -2 ), in which they 
arise from highly photoionized gas tracing the dark 
matter potentials in mildly nonlinear large-scale 
structures (for a review, see Rauch 1998). Though 
the absorbing gas is typically out of dynamical 
and thermal equilibrium, a relatively simple re- 
lation known as the Fluctuating Gunn-Peterson 
Approximation (FGPA; for the full expression see 
Croft et al. 1998) captures the essential physics 
of optically-thin Lya absorbers: 

Toc^rnr 1 , (i) 

where r is the H I optical depth, p is the local 
density of baryons (assumed to trace the dark mat- 
ter) , and Thi is the H I photoionization rate. The 
FGPA has been utilized in various forms to con- 
strain many physical parameters associated with 
high-rcdshift Lya forest absorption, such as a mea- 
surement of fif, (Rauch ct al. 1997; Weinberg et 
al. 1997) and the mass power spectrum at z ~ 3 
(Croft et al. 1998; McDonald et al. 2000a), and 
implicitly in the determination of the metallicity 
of the diffuse IGM (Rauch, Haehnelt, & Steinmetz 
1997; Dave et al. 1998) and the "equation of state" 
of the IGM (Hui & Gnedin 1997; Schaye et al. 
1999; McDonald et al. 2000b). 

At redshifts z < 1.5, the Lya transition falls in 
the ultraviolet, requiring more challenging space- 
based observation. The Hubble Space Telescope 
(HST) Quasar Absorption Line Key Project (Bah- 
call et al. 1993, 1996) obtained UV quasar spectra 
for a large sample of quasars (Jannuzi ct al. 1998) 
using the HST Faint Object Spectrograph (FOS). 
Unlike Keck/HIRES, FOS is unable to resolve in- 
dividual Lya absorption features, having a reso- 
lution of w 230 km/s. Nevertheless, this sample 
greatly expanded our understanding of the prop- 
erties and evolution of Lya forest absorbers from 
z ~ 1.7 — ► 0. One surprising result was that the 
rate of evolution of absorbers slows dramatically 
at z < 2; this was initially indicated by the large 
number of Lya lines detected in the first obser- 
vations of 3C 373 (Morris et al. 1991; Bahcall 
et al. 1991) and confirmed with a larger number 
of sight lines from the FOS Key Project (Bahcall 
et al. 1996). Cosmological hydrodynamic simula- 
tions suggested that this was a result of the dimin- 



ishing ionizing background in concert with the de- 
cline in the quasar population after z ~ 2 (Theuns, 
Leonard, & Efstathiou 1998; Dave et al. 1999, 
hereafter DHKW). FOS data also suggested that 
high-column density absorbers are evolving away 
faster than low-column density ones (Weymann et 
al. 1998), which DHKW suggested was due to 
the increase in the relative cross-section of lower- 
column absorbing gas. DHKW also showed that 
the observed correlation between absorber equiv- 
alent width and impact parameter (Chen et al. 
1998; Tripp, Lu, & Savage 1998) arises purely due 
to the clustering of matter, not due to a physical 
association of absorbers with H I in large gaseous 
halos of individual galaxies. Thus qualitatively, 
Lya absorbers at low redshift are well-described 
by current structure formation models. Still, due 
to the low resolution of FOS and the computa- 
tional difficulty of running simulations to z = 0, 
a detailed quantitative agreement of simulations 
with observations analogous to that obtained at 
high redshift has yet to be achieved. 

HST's Goddard High Resolution Spectrograph 
(GHRS) has yielded some complementary re- 
sults regarding low-z Lya absorbers. In the 
intermediate-resolution mode usually employed 
in QSO absorption line studies, the GHRS pro- 
vided an instrumental resolution of ~19 km/s, 
which is comparable to the widths of Lya lines, 
making a determination of individual absorber col- 
umn densities and 6-parameters possible, though 
not straightforward (Pcnton, Shull, & Stocke 
2000). The Far Ultraviolet Spectroscopic Explorer 
(FUSE) can be used to detect higher Lyman series 
lines at wavelengths shortward of the HST band- 
pass and thereby enable a curve-of-growth analysis 
to more accurately estimate line widths (Shull et 
al. 2000). The resulting widths were typically 
half of that obtained by direct profile fitting of 
GHRS data, which Shull ct al. (2000) interpreted 
as an indication of non-thermal components in 
Lya absorbers. 

The deployment of the Space Telescope Imag- 
ing Spectrograph (STIS) aboard HST has signifi- 
cantly improved our ability to study low-redshift 
Lya absorbers. The intermediate-resolution FUV 
echelle mode of STIS has a resolution of 7 km/s 
FWHM and enables detection of absorbers with 
^Vhi < 10 13 cm~ 2 . This resolution, similar to that 
of Keck/HIRES, allows virtually all Lya absorbers 
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to be fully resolved, facilitating a robust compar- 
ison with high-resolution hydrodynamic simula- 
tions. In this paper we determine the column 
density and linewidth distributions from STIS 
spectra of quasars PG 0953+415 (z om = 0.239) 
and H1821+643 (z cm = 0.297), and we compare 
them to carefully-constructed mock quasar spec- 
tra drawn from a cosmological hydrodynamic sim- 
ulation of a ACDM universe. We find that the 
agreement between simulations and observations 
is quite good, and lends strong support to the 
physical picture of the low-redshift intcrgalactic 
medium provided by hierarchical structure forma- 
tion models. 

DHKW suggested that the physical state of the 
gas giving rise to Lya absorption at low redshift 
is fundamentally similar to that at high redshift, 
with the exceptions that a dynamically equivalent 
patch of gas at low redshift gives rise to an ab- 
sorber with significantly lower column density (see 
their Figure 10), and collisionally ionized hot gas 
plays an increasingly important role to lower red- 
shift. With STIS data, we can now probe these 
Iow-Ahi absorbers that are predicted to be physi- 
cally similar to the high-z forest absorbers. If the 
FGPA provides a reasonable description of weak 
low-z absorbers, it may be possible to apply many 
of the same techniques used to determine physi- 
cal properties of the IGM at high redshift to STIS 
data. In this paper we present such applications, 
including an estimate of the extragalactic H I pho- 
toionization rate. Despite large systematic un- 
certainties in the modeling, our measurement is 
among the most sensitive to date. 

In §2 we describe the STIS observations of the 
two low-redshift quasars, the hydrodynamic sim- 
ulations we will use for comparison, and the con- 
struction and analysis of artificial spectra. In §3 
we present the column density and line width dis- 
tribution from STIS data alongside the distribu- 
tions from carefully-constructed artificial spectra. 
In §4 we explore the physical state of gas giving 
rise to the sorts of absorbers seen in these STIS 
spectra, and make a preliminary determination of 
the IGM temperature from the b — Nm distribu- 
tion. In §5 we constrain the average metagalac- 
tic photoionizing flux incident on these absorbers 
by matching the amplitude of the column density 
distribution to observations. We present our con- 
clusions in §6. 



2. Simulation and Observations 

2.1. STIS Quasar Spectra 

The STIS observations of PG0953+415 and 
H1821-I-643 were obtained by Tripp & Savage 
(2000) and Tripp, Savage, & Jenkins (2000) to 
study low-redshift O VI absorbers as well as the 
relationships between the various types of absorp- 
tion systems and galaxies. After the publication 
of those papers, more observations were obtained, 
and modest improvements in the data reduction 
procedures were implemented. A detailed descrip- 
tion of the observations, data reduction, and com- 
plete line lists are found in a separate paper (Tripp 
et al. 2001). Here we provide a brief summary. 

Both QSOs were observed with the medium res- 
olution FUV cchclle mode (E140M). H1821+643 
was observed with the 0'.'2 x 0'.'06 slit, while 
PG0953+415 was observed with the 0'.'2 x 0'.'2 
slit for better throughput. This STIS mode pro- 
vides a resolution of R = A/AA w 46,000 (7 km/s 
FWHM) with the Of.' 2 x 0'.'06 slit (Kimble ct al. 
1998); the FWHM is only slightly degraded with 
the OT! 2 x OT! 2 slit, but the broad wings of the line 
spread function arc more prominent (see Figure 
13.87 in the STIS Instrument Handbook, v4.1). 
This mode provides spectra extending from 1150 
— 1710 A with only a few small gaps between 
orders at A > 1630 A, redwards of the Lya forest. 

The data were reduced with the STIS Instru- 
ment Definition Team software. The individual 
spectra were flatficldcd, extracted, and wavelength 
and flux calibrated with the standard techniques. 
Then a correction for scattered light was applied 
using the method developed by the STIS Team, 
and the individual spectra were combined with 
weighting based on signal-to-noise ratio (S/N). Fi- 
nally, overlapping regions of adjacent orders were 
also coadded weighted by S/N. 

2.2. Generation of Artificial Spectra 

We employ a cosmological hydrodynamic simu- 
lation of a A-dominated cold dark matter model, 
run with PTrccSPH (Dave, Dubinski & Hernquist 
1997), having Ov = 0.6, Q C dm = 0.3527, Q b = 
0.0473, Ho — 65 kms _1 Mpc~ , inflationary spec- 
tral index n = 0.95, and rms fluctuation amplitude 
as = 0.8. We use 64 3 dark matter particles and 
64 3 gas particles in a periodic cube 11.111ft. -1 co- 
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moving Mpc on a side and a gravitational soft- 
ening length of 3.5/i -1 comoving kpc (equivalent 
Plummer softening) . Further details may be found 
in DHKW. We choose fib to match recent observa- 
tions of the deuterium abundance (Buries & Tytler 
1998), and adopt a UV background from Haardt 
& Madau (1996). In post-processing we adjusted 
the intensity of the UV background slightly in or- 
der to match the evolution of the number density 
of absorbers (dN/dz) with rest equivalent width 
W r > 0.24 A (the "matched LCDM" curve in Fig- 
ure 4 of DHKW), by assuming the optical depth 
in artificial spectra is inversely proportional to the 
photoionization rate. This simulation is otherwise 
the same as the "LCDM" run in DHKW. 

From this simulation we generate artificial spec- 
tra whose characteristics closely model those of 
PG0953+415 and H1821+643; see Dave et al. 
(1998) for a more detailed description of this pro- 
cedure. We assume zero metallicity, so all features 
in the artificial spectra are Lyman series lines. The 
flux is interpolated to the wavelength of each pixel 
in the observed spectrum, and Gaussian random 
noise is added based on the noise value at the near- 
est pixel in the observed spectrum with a similar 
flux. We convolve the spectrum with the STIS line 
spread function of the slit used to observed each 
quasar, interpolated to the pixels' wavelength. We 
fit a continuum to both observed and artificial 
spectra using a median filter in a sliding 100-pixel 
(« 300 km/s) window. We adjust the median fil- 
ter so that in each set of artificial spectra, there is 
on average no net flux subtracted by the contin- 
uum fitting (i.e. we recover the true continuum, 
on average). We generate 20 artificial spectra for 
each quasar, for a total of 40. 

We fit Voigt profiles to the observed and arti- 
ficial spectra using AutoVP (Dave et al. 1997). 
We choose a 4.5cr detection threshold for identify- 
ing Lya absorbers, using the procedure described 
in Lanzetta, Turnshek, & Wolfe (1987) with a 
20 km/s detection window. We only consider lines 
between the Lya and Ly/3 emission peaks so as to 
reduce confusion with higher-order Lyman series 
and metal lines, and do not include any absorbers 
within 5000 km/s of the Lya emission peak. Spec- 
tral regions containing Galactic ISM lines or ex- 
tragalactic metal lines identified in the data are 
discarded, and any lines identified in the artifi- 
cial spectra in those wavelength regions are dis- 



carded as well, ensuring that each artificial spec- 
trum has the same redshift path length as its cor- 
responding quasar spectrum. Note that consider- 
able care must be exercised in the identification of 
the Milky Way features since the spectra have suf- 
ficient sensitivity to show weak ISM lines such as 
the C I and Ni II resonance transitions; similarly, 
several extragalactic metal lines are present in the 
wavelength range of interest (see Tripp, Savage, & 
Jenkins 2000; Tripp et al. 2001). 

An 80A sample of H1821+643 is shown in Fig- 
ure 1, at three stages. The first (top panel) shows 
the unnormalized observed spectrum, with the 
continuum fit as the thick solid line, and the la- 
noise vector at the bottom. The middle panel 
shows the normalized spectrum with Voigt profile 
identifications indicated by the tick marks (some 
of the marked lines are Milky Way or extragalac- 
tic metal lines which are discarded). The bot- 
tom panel shows the corresponding segment from 
an artificial H1821+643 spectrum, with the same 
analysis procedure having been applied. 

While there may be biases associated with our 
particular continuum fitting, line identification, 
and Voigt profile fitting algorithms, the applica- 
tion of the identical routines to both simulations 
and observations facilitates a fair comparison be- 
tween the two. We do not expect these biases to 
be significant since we arc fully resolving virtu- 
ally all Lya absorbers, and continuum fitting is 
relatively straightforward since the low-z forest is 
sparse. Any systematic biases that do arise will, 
in most cases, be reflected equally in the statistics 
of both the simulations and observations. 

Our ll.lll/i _1 Mpc simulation volume is small 
compared to the scale of linearity at z = 0, hence 
our derived statistics are subject to cosmic vari- 
ance. We do not attempt to estimate this, since 
to do so would require running a suite of simula- 
tions that is beyond our current capabilities. We 
note that our randomly chosen volume is typical, 
containing no anomalously large mass concentra- 
tions or voids. As we will show, most Lya ab- 
sorbers arise in moderate overdensity regions that 
are well-sampled even within our small volume, 
though the contribution from coupling to waves 
larger than our boxsize is not included. We leave 
a systematic study of these effects for future work. 
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3. Absorber Statistics 




Fig. 1. — An 80 A segment from the Lya forest of 
H1821+643. Top panel shows the unnormalized 
spectrum with continuum fit. Middle panel shows 
the normalized spectrum with line identifications 
by AutoVP. Note that some of the marked lines 
in this panel are not Lya lines, e.g., the Si III 
1206.5A lines at z = 0.225 (see Figure 2 in Tripp, 
Savage, & Jenkins 2000) and the Milky Way Ni II 
line at 1454. 8A. Bottom panel shows an artificial 
HI 82 1+643 spectrum mimicking its resolution and 
noise characteristics, continuum fit and analyzed 
with the same routines applied to the H1821+643 
spectrum. The ler noise vector is shown near the 
bottom of each panel. 



The distributions of column densities and 
Dopplcr widths b obtained by fitting Voigt pro- 
files to spectral features have historically been 
used to characterize the statistics of Lya forest 
absorbers. Though these statistics reflect the old 
paradigm that Lya absorbers arise from isolated 
thermally-broadened clouds, an assumption that 
is likely to be incorrect in detail (e.g. Hernquist 
et al. 1996; Outram, Carswell & Theuns 2000; 
Theuns, Schaye, & Hachnelt 2000), such statistics 
still represent a fair statistical characterization of 
the absorber population, more so at low redshift 
than at high redshift because line blending is less 
severe and individual absorbers can more easily 
be identified. 

At high redshift, statistics based on the distri- 
bution of flux in the Lya forest have gained fa- 
vor (e.g. McDonald et al. 2000a; Machecek et al. 
2000) due to their simplicity, ease of implemen- 
tation, and robustness of comparison with mod- 
els. At low redshift, these statistics are somewhat 
more difficult to apply reliably, because they re- 
quire that the opacity in the forest be predom- 
inantly due to Lya absorption, whereas at low- z 
the relative contribution from metal and ISM lines 
is greatly increased. Thus we will leave flux statis- 
tics for future work, and focus here on Voigt profile 
statistics. 

We identify 33 Lya absorption lines in the Lya 
region of PG0953+415, and 56 in H1821+643, 
using AutoVP. There arc more absorbers in 
H1821+643 because the average signal-to-noise 
ratio per w 3 km/s pixel in this wavelength 
range is higher, S/N= 13.7 compared to 8.1 
for PG0953+415, so more weak absorbers are 
indentified. The total redshift path length is 
Az = 0.165 for PG0953+415 and Az = 0.170 
for HI 82 1+643, not including discarded regions 
that cover Az = 0.008 and Az = 0.011, respec- 
tively. This path length is for absorbers above 
our completeness limit, which we will show is 
A^hi ^ 10 13 cm -2 ; for weaker absorbers, the effec- 
tive path length is smaller (see, for example, the 
discussion in Penton, Shull, & Stocke 2000), but 
by construction is still identical in the data and 
the artificial spectra. We will only quote statis- 
tics for absorbers above our completeness limit. 
The mean redshift of our entire absorber sample 
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is z = 0.17. 



3.1. Column Density Distribution 



The column density distribution is denned as 
the number of lines per unit redshift z per unit 
column density Ajji, and is generally parameter- 
ized as a power law in Nm, 



f(Nm) 



d 2 N 
dzdNm * HI 







(2) 



At z > 2.5, (3 w 1.5 (Hu ct al. 1995; Lu ct al. 
1996; Kirkman & Tytler 1997; Kim et al. 1997). 
Recent work with VLT/UVES suggests a similar 
slope ((3 « 1.4) down to z ~ 1.5 (Kim, Cris- 
tiani & D'Odorico 2000). This is predicted to 
steepen significantly at low redshift (Zhang, An- 
ninos & Norman 1995; Theuns, Leonard, & Ef- 
stathiou 1998, DHKW). GHRS observations by 
Penton, Shull, & Stockc (2000) indicate a some- 
what steeper slope of (3 = 1.8, but that was ob- 
tained by assuming a particular b for all absorbers 
(varied between 20 and 30 km/s) rather than by 
determining b and TVhi independently from profile 
fitting. An alternative statement of this trend is 
that high-A^Hi absorbers are predicted to evolve 
away faster than Iow-A^hi ones. The analysis in 
§5.2 of Penton, Shull, & Stocke (2000) suggests no 
significant difference between the evolution of ab- 
sorbers with 10 131 < N m < 10 14 cuT 2 and N m > 
10 14 cm~ 2 , contradicting FOS results (Weymann 
ct al. 1998) as well as model predictions. How- 
ever, Figure 18 of Penton, Shull, & Stockc (2000) 
suggests that between 10 13 < A H i < 10 14 cm~ 2 , 
stronger absorbers evolve away significantly faster 
than weaker ones from z = 3 — > 0. 

Figure 2 shows the column density distribution 
of identified absorbers in the two STIS spectra 
(combined) along with the corresponding results 
from 40 artificial spectra. The simulations and ob- 
servations show very good agreement, in both the 
slope and the amplitude of f(Nui)- The best- fit 
slopes for N m > 10 13 c m - 2 are f3 ohs = 2.04 ± 0.23 
and (3 s i m — 2.15 ± 0.04. Thus we find that the the 
column density distribution has steepened consid- 
erably since z ~ 3, at least for absorbers with 
10 13 < Ahi < 10 14 cm~ 2 , in general agreement 
with Figure 18 of Penton, Shull, & Stocke (2000). 
We do find a slope that is somewhat steeper (by 
~ la) than that found from the GHRS sample 
(Penton, Shull, & Stocke 2000). Note that we 
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Fig. 2. — Column density distributions from ob- 
servations of PG0953+415 and H1821+643 (solid 
line) and a sample of 20 artificial spectra for each 
quasar (dashed line) . Dotted line shows the result 
from artificial spectra with S/N increased by x5, 
indicating that the absorber population is com- 
plete down to at least TVhi ~ 10 12 - 9 cm~ 2 . Note 
that the ionizing background used to generate the 
artificial spectra was normalized to the evolution 
of absorbers in the Key Project FOS sample (see 
DHKW), and is not the the normalization that 
best fits this STIS sample; we will quantify this 
further in §5. 
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only include ./Vhi > 10 12 ' 9 cm~ 2 bins in this fit; 
had we included the bin at Nm = 10 12 ' 6 cm -2 , the 
slopes would reduce to « 1.7. This bin is only 
slightly incomplete (as we discuss below), but it 
illustrates that completeness must be carefully as- 
sessed before an accurate determination of /3 can 
be made. 

The dotted line in Figure 2 shows the column 
density distribution for a set of artificial spectra 
where the S/N ratio has been increased by a fac- 
tor of 5, allowing weaker absorbers to be identi- 
fied. No significant differences are seen for Ahi > 
10 12 ' 9 cm~ 2 . This tests the completeness of the 
observations (and simulations) at low JVhi, show- 
ing that the absorber sample in PG0953+415 and 
H1821+643 is complete down to N m « 10 13 cnT 2 , 
and > 70% complete to N m « 10 12 - 6 cm" 2 . The 
similar shape of the turnover below this value in- 
dicates that the noise level in the artificial spectra 
are faithfully mimicking that in the STIS data. 

The agreement in amplitude of /(-/Vhi) de- 
pends on the chosen value of the H I photoion- 
ization rate used to generate the artificial spec- 
tra, as we will explore further in §5. The best- 
fit amplitudes for the simulated and observed ab- 
sorbers are / sim (10 13 cm" 2 ) = -10.70 ± 0.02 and 
/ Obs (10 13 cm- 2 ) = -10.87 ±0.12 respectively. Re- 
call that the ionizing background was set to obtain 
agreement with dN/dz of absorbers with W r > 
0.24A (N m > 10 14 cm- 2 ) measured by the Key 
Project's sample of FOS quasar spectra (see Dave 
et al. 1999), which probes a different column den- 
sity regime than these STIS spectra. The general 
agreement between these independent data sets 
lends further support to the overall model. 

3.2. Distribution of 6-parameters 

The 6-parametcr distribution has historically 
provided a significant challenge for theories of the 
Lya forest. At high redshift, absorbers are pre- 
dominantly broadened by bulk flow (e.g. Hern- 
quist et al. 1996; Outram, Carswell & Theuns 
2000), but thermal broadening is not negligible 
(Theuns, Schaye, & Haehnelt 2000). Hence the 
6-parameter distribution can probe the tempera- 
ture of Lya forest gas (Schaye et al. 1999; Mc- 
Donald et al. 2000b; Ricotti, Gnedin & Shull 
2000). The derived temperatures seem to sug- 
gest that intergalactic gas at z ~ 3 is heated 
by some other process in addition to equilibrium 



H I photoionization, such as helium reionization 
(Haehnelt & Steinmetz 1998). Though some ob- 
servational constraints are beginning to emerge 
(Heap et al. 2000; Songaila & Cowie 1996; Giroux 
& Shull 1997), theories for when and how reion- 
ization occurs are quite uncertain, making an a 
priori prediction of 6-parameters difficult. Fur- 
thermore, it has been shown that a robust predic- 
tion of the 6-parameter distribution for compari- 
son with Keck/HIRES data requires exceptionally 
high numerical resolution (Theuns et al. 1998; 
Bryan et al. 1999), higher for instance than in 
the simulation analyzed here. 

At low redshift, it is unclear whether simula- 
tion predictions of 6-parameters are more or less 
robust. Certainly reionization is far enough in the 
past that the IGM is expected to have returned 
to an equilibrium ionization condition. However, 
the character of the low-z IGM has changed dra- 
matically. In particular, though a substantial frac- 
tion of the baryons remains in cool, photoionized 
gas, a significant component of intergalactic gas 
now resides at temperatures exceeding that from 
photoionization, due to shock heating from infall 
onto large-scale structure (DHKW; Cen & Os- 
triker 1999; Dave et al. 2001, hereafter D01). Fur- 
thermore, it is possible that non-thermal processes 
such as supernova heating also heat the IGM (Wu, 
Fabian, & Nulsen 2000, though see D01 and Croft 
et al. 2000). The simulation used here includes 
these physical processes (it is simulation "D2" in 
D01), but they are more complex than just pho- 
toionization heating and adiabatic cooling, involv- 
ing highly uncertain theoretical issues such as es- 
cape fraction of feedback energy and the contri- 
bution of small-scale shocks to IGM heating (see 
D01 for discussion). Furthermore, the baryonic 
mass fraction in various intergalactic components 
depends on and cosmology (DHKW) in a way 
that is difficult to predict analytically, as we will 
discuss in §5. 

With these caveats in mind, we present in Fig- 
ure 3 a comparison of the 6-parameter distribu- 
tion in our two STIS quasar spectra and the 
artificial spectra. We only show the distribu- 
tion for absorbers above our completeness limit of 
Ahi > 10 13 cm~ 2 , but including absorbers down 
to Nm = 10 12,6 cm~ 2 makes little difference in 
the results. Remarkably, despite all the poten- 
tial pitfalls, the agreement between simulations 



7 



0.1 



0.08 r 



§0.06 
&).04 



0.02 




i — 1 — 1 — 1 — i — 1 — 1 — 1 — r 
-, Observed — 

ACDM - 

Thermal Only ■■ 



20 40 60 

b (km/s) 



Fig. 3. — ^-parameter distribution from obser- 
vations (solid line) and artificial spectra (dashed 
line). The dotted line shows the distribution of 
thermal widths for the absorbers in the artifi- 
cial spectra, discussed in §4.3. Only lines with 
N m > 10 13 cm~ 2 are included. 



and observations is quite good. The median b- 
parameter is observed to be 22 km/s, while the 
simulations predict 21.5 km/s. The average b 
is observed to be 25 km/s, and predicted to be 
23 km/s. A K-S test shows no statistically signif- 
icant difference between the b distributions from 
the observations and simulations. These quoted 
values include the contribution from the STIS line 
spread function; correcting for this lowers these 
values by ~ 1 km/s, yielding an intrinsic median 
^-parameter of w 21 km/s. The dotted line shows 
the distribution of thermal widths for the artificial 
sample, which we will discuss in §4.3. 

The median b value is lower than the 27 — 
36 km/s typically found at z ~ 2.5 — 4 (Kim et al. 
1997), though a higher S/N spectrum finds some- 
what lower b values (Kirkman & Tytler 1997). 
This is qualitatively consistent with the predic- 
tion that the unshocked IGM is cooler at low- 2 
due to the lowered intensity of the photoionizing 
background, but could also be a result of lowered 
bulk flow broadening due to differential Hubble ex- 
pansion across the absorbing structures. We will 
explore the relative importance of these effects in 
§4.3. 

Our 6-parameters do not follow an extrapo- 
lation of the trend seen by Kim et al. (1997), 
that b generally increases with decreasing redshift 
for 2.1 < z < 4. They are also significantly 
smaller than ^-parameters derived from profile fit- 
ting GHRS data (b = 38 ± 16 km/s, Penton, Shull, 
& Stocke 2000). On the other hand, we are in 
agreement with a curve-of-growth analysis using 
FUSE and GHRS data by Shull et al. (2000), who 
find b = 31±7 km/s and 6 mo dian = 28 km/s. As we 
will discuss further in §4.3, larger column density 
absorbers are wider, and Figure 6 shows that at 
the median JVhi = 10 14,2 cm~ 2 of the Shull et al. 
(2000) sample, our median 6-parameter is ps 30. 
This agrees quite well with their value of 28 km/s, 
given the small samples involved. Our STIS data 
probes weaker lines than the FUSE sample, and 
hence results in a smaller median b. 

Shull et al. (2000) find that, for partially sat- 
urated lines, direct profile fitting of GHRS data 
overstimates the ^-parameter by a factor of two 
on average, and thereby calls into question the va- 
lidity of profile fitting. It is true that results from 
profile-fitting a saturated Lya line can have large 
uncertainties. However, the vast majority of our 
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lines are not saturated (Ahi < 10 14 cm~ 2 ), so our 
statistics, particularly our median b value, should 
not be too sensitive to this effect. Furthermore, 
as evidenced by our agreement with the curve- 
of-growth study of Shull et al. (2000), we are 
at most only slightly overestimating 6-parameters. 
The 7 km/s resolution of our data insures that we 
can adequately sample individual Lya absorbers 
profiles, since the smallest Lya absorber widths 
are significantly larger. Note that we do find a 
number of metal and ISM lines that have widths 
comparable to the resolution (Tripp, Savage, & 
Jenkins 2000; Tripp et al. 2001), indicating that 
if such Lya lines existed, they would have been 
detected. 

The agreement of the predicted and observed 
6-parameter distributions is actually better than 
that obtained in a similar comparison at high 
redshift (Dave et al. 1997). This surprising re- 
sult may be because STIS observations do not 
yet have high enough S/N to show discrepancies 
with the simulations. Or it may be that since 
low- 2 forest requires less deblending of absorption 
features, the 6-parameter distribution is more ro- 
bustly predicted. Alternatively, we may have sim- 
ply been fortunate with this simulation, and fu- 
ture higher-resolution, larger-volume simulations 
will show discrepancies. For now we are encour- 
aged by the good agreement, and will proceed on 
the assumption that it is not a spurious coinci- 
dence, but rather reflects the basic validity of the 
underlying physical model. 

4. Physical Properties of Lya Absorbers 

The good agreement between the observed and 
simulated absorber statistics suggests that the 
simulations are accurately describing the origin 
of weak Lya absorbers at low redshift. In this 
section, we investigate the physical conditions in 
the simulated IGM that give rise to weak Lya ab- 
sorbers. These studies are similar in spirit to those 
in DHKW, but with a focus on the weak absorber 
population seen in these STIS spectra. We also in- 
vestigate the relationship between the 6-parametcr 
distribution and the temperature of the IGM. 

4.1. The Temperature-Density Relation 

At high redshift, there is a tight relationship 
between density p and temperature T in the gas 



producing the Lya forest, arising from the bal- 
ance between photoionization heating and adia- 
batic cooling due to Hubble expansion (e.g. Hui 
& Gnedin 1997). At low redshift, shocks have 
heated a significant fraction of even diffuse inter- 
galactic baryons above the equilibrium photoion- 
ization temperatures, diluting this tight relation 
(e.g. Cen & Ostriker 1999, DHKW). However, 
this shock-heated gas tends to have a low ioniza- 
tion fraction, and hence is less likely to give rise 
to Lya absorbers, while a substantial amount of 
baryonic matter remains in the cool, diffuse IGM 
(see Figure 1 in D01). Thus Lya absorbing gas 
may still have a reasonably tight p — T relation. 

Figure 4 shows the temperature T vs. density 
p/p of all absorbers identified in our 40 artificial 
spectra. The density and temperature of an ab- 
sorber are taken to be those at the location of the 
absorber's maximum optical depth, and the den- 
sity is shown in units of the mean baryonic density 
at the absorber redshift. A significant fraction of 
the absorbers lie along a relatively tight relation in 
p—T space arising from photoionization, but there 
are also many shock-heated absorbers at higher T. 
The photoionized absorbers follow the fitting for- 
mula 

T« 5000 (J J K, (3) 

shown as the dashed line, whose slope is similar to 
that seen at high redshift (Croft et al. 1998). For 
comparison, we show contours encompassing 50% 
and 90% of all gas in our simulation. In the den- 
sity range of STIS Lya absorbers, the absorbers' 
range of temperatures is significantly smaller than 
that of all baryons. Evidently, Lya absorbers pref- 
erentially (though not uniquely) trace out gas in 
an equilibrium photoionized state. 

4.2. Column Density vs. Gas Density 

The FGPA indicates a power-law relationship 
between the optical depth and the underlying gas 
density. At low redshift, this relationship becomes 
contaminated by hot gas that lies off the p — T 
relation as seen in Figure 4, resulting in absorbers 
with more widely varying ionization fractions, and 
hence more widely varying underlying densities. 

Figure 5 show plots of column density Ahi vs. 
density of underlying gas for the absorbers iden- 
tified in the 40 artificial spectra. As suggested 
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Fig. 4. — Temperature vs. density for absorbers 
from our artificial spectra. The contours enclose 
50% and 90% of all gas particles in our simulation. 
These contours indicate that while an important 
shock-heated population is present, the Lya ab- 
sorbers preferentially sample gas in the photoion- 
ized regime. The dashed line shows a fit to the 
photoionized absorbers, given in equation (3). 



Fig. 5. — Column density JVhi vs. density in units 
of the mean p/p, for absorbers in our artificial 
spectra. Left and right panels show results from 
H1821-like and PG0953-like artificial spectra, re- 
spectively. The line shows the relation given in 
equation (4). 
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by Figure 10 of DHKW, there is substantially 
greater scatter at low redshift as compared to high 
redshift. The line shows a similar fit to that in 
DHKW, 



P 
P 



12 



HI 



10 14 cm- 2 



0.7 



10 



-0.42 



(4) 



corrected for the difference in Tri used (discussed 
further in §5). This is a reasonable fit to the ar- 
tificial ST1S absorbers, though there is a hint of 
greater scatter than in the z = absorber pop- 
ulation of DHKW. This may be because DHKW 
used S/N=30 per pixel, whereas our STIS spectra 
have S/N~ 6—15 per pixel, resulting in larger er- 
rors in profile fitting to obtain Nhi, especially at 
lower column densities. In Figure 5, the left and 
right panels show the results from H1821-like and 
PG0953-like artificial spectra, respectively. The 
most deviant outliers in density (at low column 
densities) arise in the lower-S/N PG0953-like spec- 
tra, indicating that S/N issues may be important. 

Given the large scatter, it is difficult to unam- 
bigously assign an underlying density to an ab- 
sorber of a given Nui, though a trend clearly ex- 
ists. While higher quality data may tighten this 
relation somewhat, as of now it appears the FGPA 
is not nearly as tight a relation at low redshift as 
compared to high-z. Thus techniques such as re- 
covery of the mass power spectrum (Croft et al. 
1998) will face greater challenges in this regime. 

4.3. The Temperature of the IGM 

Since higher column density systems arise from 
higher density gas, and higher density gas has 
higher temperatures, one would expect that b 
should be correlated with Am- In particular, 
the smallest 6 parameters at a given Nhi should 
mainly reflect the temperature of the underlying 
absorbers, as thermal broadening should dominate 
in the narrowest lines. This fact has been used 
at high redshift to constrain the temperature and 
"equation of state" (i.e. the p — T relation, Hui 
& Gnedin 1997) of the IGM (Schaye et al. 1999; 
McDonald et al. 2000a), and a similar argument 
has been extended to low redshift as well (Ricotti, 
Gnedin & Shull 2000). It is conceivable that re- 
covering the IGM temperature is actually easier 
at low-z, because at high redshift, narrow lines 
can arise from Voigt profile fitting the asymmet- 
ric wings of larger lines, and may be confused for 



true narrow absorbers. This occurs much less fre- 
quently in the sparser low redshift forest. 

Figure 6 shows a plot of Nui vs. b for the ab- 
sorbers identified in PG0953+415 and H1821+643 
(left panel) and absorbers in our artificial spectra 
(right panel). The dashed curve shows a running 
median of b parameters with iVm, with variance 
in each bin. There is a clear trend for larger b- 
parameters at larger column densities, as would 
be expected by the above argument. 

We can quantify this relation as follows. In 
§4.2 we confirmed that equation (4) is a reason- 
able fit to STIS absorbers (though with large scat- 
ter). Equation (3) gives our p — T fit to pho- 
toionized absorbers, which are the coolest at any 
given density. Combining these relations and us- 
ing b = y // 2kBT/m p for pure thermal broadening, 
we obtain an estimate for the minimum b param- 
eter as a function of column density and redshift: 



^thermal 



19 



HI 



0.21 



10 



-0.12x 



km/s. 



(5) 

This relation is shown as the solid line in the two 
panels of Figure 6, for z = 0.17 (the mean red- 
shift of the observed sample). It is a good ap- 
proximation to the lower envelope of absorber b- 



parameters for A H i > 10 



13.3, 



in our artificial 



spectra (right panel). Below this column density, 
presumably blending is more frequent, and narrow 
lines from line wings dilute this relationship. Also, 
noise can cause some weak lines to appear nar- 
rower than they really are, since weaker lines have 
larger profile fitting uncertainties. As seen in the 
left panel, equation (5) is consistent with a lower 
envelope for all but the weakest STIS absorbers as 
well, though a larger sample of absorbers will be 
required before such a relation may be recovered 
directly from the observations. 

Comparing the solid line (pure thermal broad- 
ening) with the dashed line (median b) in Figure 6 
suggests that thermal broadening is responsible for 
a substantial portion of the absorber widths, with 
^thermal ~ 0.76 for a typical absorber. This frac- 
tion remains fairly constant with Am- This is also 
seen in Figure 3, where the dotted line shows the 
histogram of thermal widths for Am > 10 13 cm~ 2 
absorbers; its median value is ~ 15 km/s. The 
contribution from thermal broadening is there- 
fore greater than at high redshift, where bulk flow 
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Fig. 6. — Nhi vs. 6-parameter for STIS-observed 
absorbers (left panel) and absorbers in 40 artificial 
spectra (right panel) . The thick dashed line shows 
the median b in four bins of Nm, with variance 
computed within each bin. The solid line shows 
thermal broadening as a function of iVm, given by 
equation (5). 
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broadening dominates (e.g. Hernquist et al. 1996; 
Schaye et al. 1999). 

Equation (5) was obtained from fits to the par- 
ticular simulation used here; it is not true in gen- 
eral. In order to measure the temperature of the 
IGM unambiguously, it is important to test equa- 
tions (3) and (4) over a wider range of parameter 
space. Nevertheless, the good agreement with the 
6-parameter distribution, as well as the significant 
contribution of thermal broadening to this distri- 
bution, suggests that these simulations are accu- 
rately characterizing the "equation of state" of the 
low-2 IGM. Our measurement is given by equa- 
tion (3), indicating a temperature of <~ 5000 K at 
the mean density. Increasing the temperature by 
more than a factor of two would produce thermal 
6-parameters alone that are larger than observed, 
so we place an upper limit at T(p) < 10 4 K. Since 
this fit is only for our current simulation, it is sub- 
ject to the uncertainties in modeling, but at least 
at high redshift measurements of T are not very 
sensitive to such uncertainties (Theuns, Schaye, & 
Haehnelt 2000). Note that we are measuring here 
the temperature of equilibrium photoionized Lya 
forest absorbers only, and this does not exclude 
the possibility that the IGM contains a compo- 
nent of hotter gas that may give rise to, say, O VI 
absorption (Tripp, Savage, & Jenkins 2000). 

A similar technique has already been used 
to determine the temperature of the IGM from 
GHRS data, obtaining a value of T(p) = 4700 K 
for J = 10 _23 ergs _1 cm~ 2 sr _1 Hz -1 (Ricotti, 
Gnedin & Shull 2000), in agreement with our re- 
sult. With more STIS data and better simulations, 
this type of analysis should be able to robustly 
measure the temperature of the low-redshift IGM 
in the near future. 

5. Preliminary Constraints on Thi 

As stated in §3, the agreement between the ob- 
served and simulated column density distribution 
amplitude is in principle a test of the ionization 
state of the absorbing gas. Because the absorb- 
ing gas for 7V H i < 10 14 cm~ 2 is optically thin, 
Nyu oc T H i \ so changing the ionization 
fraction of Lya forest absorbers will shift the simu- 
lated f(N}n) in the horizontal direction, and alter 
the agreement with observations. 

It is straightforward to determine the shift in 
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TVhi that produces the best agreement with obser- 
vations, then translate that into a value for Thi- 
This yields T H i = l(r 13 - 3±0 - 12 s" 1 at an average 
rcdshift of z = 0.17, up by « 0.2 dex from the 
value of Thi originally used to construct the ar- 
tificial spectra. Figure 7 shows this value of Fhi 
(filled circle) in relation to other predicted and ob- 
served values. The error bars on this value are dis- 
cussed below. The value derived here is close to 
the "matched LCDM" curve from DHKW (solid 
line), since it was this J v that produced the good 
agreement with the observed f(Nm) seen in Fig- 
ure 2. This value is even closer to that predicted 
from the evolution of the quasar population by 
Haardt & Madau (1996, dashed line) and three- 
fourths of a similar prediction by Fardal, Giroux, 
& Shull (1998, dotted line) at z = 0.17. A recent 
calculation accounting for emission from QSOs, 
Seyferts and starbursts (Shull et al. 1999) is shown 
as the diamond at z = 0, and is quite consistent 
with our measurement. 

Figure 7 also shows various measurements that 
constrain Thi at low redshift (a more complete dis- 
cussion may be found in Shull et al. 1999). The 
open circle at z ~ 0.5 comes from the proximity 
effect (Kulkarni & Fall 1993), the open square 
comes from modeling the radial truncation of H I 
in the isolated spiral NGC 3198 (Maloney 1993; 
Dove & Shull 1994), and the two downward tri- 
angles show upper limits from Ha emission due 
to photoevaporation off nearby H I clouds (Vogel 
et al. 1995; Donahue, Aldcring & Stocke 1995). 
An actual measurement using this last technique 
is in preparation (R. Weymann, private commu- 
nication). Tumlinson et al. (1999) used metal 
absorption line ratios in the halo gas of NGC 3067 
to model the required photoionizing background, 
and their lower limit is shown as the upward trian- 
gle. Note that all these methods assume that the 
ionizing radiation is mctagalactic and not local, 
despite measuring the field close to an observable 
object. Conversely, weak Lya absorbers are typi- 
cally far from galaxies (Stocke et al. 1995; Tripp, 
Lu, & Savage 1998). Additionally, our statistical 
errors are smaller than those of any of the above 
measurements, but as we now discuss, our model- 
ing uncertainties dominate. 

The FGPA (eq. 1) indicates that r depends on 
the density of the absorbing gas in addition to r H i- 
The exponent of this dependence is approximately 
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Fig. 7. — Our measurement of Thi (filled circle) 
from the column density distribution of the Lya 
forest. The errors are dominated by modeling un- 
certainties; the horizontal errors show the redshift 
range of absorbers in our STIS spectra. Predic- 
tions for the evolution of Thi are shown as dashed 
(Haardt & Madau 1996) and dotted (Fardal, 
Giroux, & Shull 1998) lines. Solid line shows the 
"matched LCDM" rate that reproduces the evo- 
lution of dN/dz for W r > 0.241 FOS absorbers 
(DHKW). This Thi was the one used to generate 
our artificial STIS spectra. Various measurements 
are indicated: proximity effect (open circle, Kulka- 
rni & Fall 1993); H I truncation in spiral galaxy 
(open square, Maloney 1993; Dove & Shull 1994); 
photoevaporation off H I cloud (downward trian- 
gles, upper limits, Vogel et al. 1995; Donahue, 
Aldering & Stocke 1995); modeling of metal line 
ratios (upward triangle, lower limit, Tumlinson et 
al. 1999). Some points near z = have been offset 
slightly for clarity. 
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1.6 at high redshift. In the previous section we 
showed that low redshift absorbers are consistent 
with this relation, though the spread is signifi- 
cantly greater. What is less certain is how the 
density of the absorbing gas varies as one changes 
fib, cosmology, and numerical parameters in the 
simulations. 

At high redshift, nearly all the baryons are in 
the Lya forest, so the FGPA can be used to di- 
rectly probe O;, (Rauch et al. 1997; Weinberg et 
al. 1997). At low redshift, the fraction of gas giv- 
ing rise to Lya absorbers is significantly smaller 
(~ 20%; Penton, Shull, & Stocke 2000), because 
many of the baryons have moved into other phases 
(DHKW, Cen & Ostriker 1999, D01). Though 
modeling uncertainties are large (as we discuss be- 
low), the cosmological simulations generally pre- 
dict that at the present epoch roughly 1/3 of the 
baryons are in diffuse photoionized gas, 1/3 are 
in the shock- heated IGM, and 1/3 are in virial- 
ized systems such as galaxies and clusters (see D01 
and references therein). Studies of low-z Lya ab- 
sorbers suggest that the diffuse photoionized IGM 
contains a considerable fraction of the baryons 
(Shull, Stocke & Penton 1996; Penton, Shull, & 
Stocke 2000). The shock-heated IGM is more 
challenging to detect, but detections of O VI at 
low redshift suggest that the predicted hot gas 
is indeed present (e.g. Tripp, Savage, & Jenkins 
2000), though since O VI absorption can occur in 
photoionized gas and Lya lines can originate in 
hot gas (see, e.g., §4 in Tripp & Savage 2000), 
there may be some double-counting in estimates 
of baryon reservoirs traced by various QSO ab- 
sorbers. 

Exactly how much baryonic matter remains in 
gas traced by Lya absorbers is uncertain. One es- 
timate may be obtained if we consider that low-z 
absorbers preferentially arise in the so-called "dif- 
fuse" (i.e. photoionized) phase of DHKW and 
D01, as suggested by Figure 4. In D01, the baryon 
fraction in this phase varies by as much as a fac- 
tor of two (~ 20 — 40%) among their simulations 
having various volumes, resolutions, and hydrody- 
namic algorithms. Furthermore, those simulations 
all were based on ACDM universes with similar f^. 
While this world model has gained much support 
recently (Bahcall et al. 1999), DHKW's analy- 
sis of several different cosmologies also indicates 
a possible variation of a factor of <~ 1.5 in the 



baryon fraction in this phase due to the differing 
rates of structure formation. Additionally, there 
is some uncertainty in fib itself, and uncertainties 
in Ho come into play as well since the observed 
quantity is usually the combination f^/i 2 . Due 
to these variations, combined with r oc p 1 - 6 /Tni, 
our best estimate for a la systematic uncertainty 
in the determination of Thi is roughly a factor of 
~ 5. Clearly this dominates over the statistical 
error quoted above. 

Thus we will claim a preliminary measurement 
of r ffl = 10- 13 - 3±0 - 7 S - 1 from the Lya forest, at 
z = 0.17. This value is shown as the filled cir- 
cle in Figure 7; the horizontal error bar indi- 
cates the redshift range covered by our STIS spec- 
tra. For the spectral shape proposed by Fardal et 
al. (similar to that of Haardt & Madau and ap- 
proximately v -1 - 8 ), this translates to J„ = 1.9 x 
10~ 23 erg s -1 cm~ 2 sr _1 Hz -1 at 912A. Our prelim- 
inary result would clearly benefit from improved 
simulations to better quantify the uncertainties in- 
volved. Unfortunately, analytic or semi-analytic 
methods such as that of Hui, Gncdin, & Zhang 
(1997) are not applicable to this problem, since the 
physics of the diffuse IGM at low-z involves a com- 
plex dynamical interplay of non-equilibrium struc- 
tures. A proper quantification would involve per- 
forming a suite of large-scale hydrodynamic sim- 
ulations exploring the relevant parameters, some- 
thing that is beyond our current computational ca- 
pability. Still, despite the large and poorly quan- 
tified uncertainties in the value derived here, ours 
is among the best observational constraints on the 
metagalactic H I photoionization rate to date, and 
offers hope that a more precise determination may 
be forthcoming soon. 

6. Conclusions 

We have examined the intergalactic Lya ab- 
sorber population in STIS spectra of PG0953+415 
and H1821-I-643, and compared it to that in arti- 
ficial spectra drawn from a cosmological hydrody- 
namic simulation of a ACDM universe. The good 
spectral resolution (7 km/s) and relatively high 
S/N of these STIS spectra yield an unprecedented 
view into the weak absorber population at low red- 
shift. We find: 

• The column density distributions agree quite 
well in slope and amplitude. The measured 
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slope is (3 — 2.04±0.23, considerably steeper 
than that seen at high rcdshift (fi w 1.5), in- 
dicating that strong absorption line systems 
have evolved away faster than weak ones. 

• The 6-parameter (linewidth) distributions 
agree quite well, with an intrinsic median 
6-parameter of ~ 21 km/s. This value is 
lower than that seen at high rcdshift, sug- 
gesting that the IGM has cooled and/or bulk 
flow broadening has lessened. The median 
6-parameter increases with TVhi, indicating 
that denser IGM gas is hotter. From a com- 
parison of 6-parameters with the expected 
thermal linewidth, thermal broadening has 
an increased contribution to linewidths at 

10W-Z, With ^thermal ~ 0.76. 

• Lya absorbers preferentially arise in gas hav- 
ing temperatures close to that expected from 
photoionization, with T oc p o e , though there 
is substantial scatter to higher temperatures 
due to absorbers arising in the shock-heated 
IGM. 

• The density-column density relation, a.k.a. 
the fluctuating Gunn-Peterson approxima- 
tion, is considerably less tight at low red- 
shift than at high-z. A trend still exists, 
consistent with that seen at high redshift, 
but shows greater scatter. 

• The minimum 6-parameter as a function 
of column density faithfully yields the 
temperature-density relation in the IGM. 
The simulation analyzed here suggests T ~ 
5000 K, with an upper limit of 10 4 K, for 
Lya absorbing gas at the cosmic mean den- 
sity. 

• The amplitude of the column density dis- 
tribution may be used to place constraints 
on Thi, the metagalactic H I photoioniza- 
tion rate at low-z. We find, at z « 0.17, 
r H i ~ io- 13 - 3±0 - 7 s _1 , with systematic un- 
certainties in modeling being the dominant 
source of error. 

In summary, weak Lya absorbers in the low- 
redshift IGM seen with STIS appear to be physi- 
cally similar to high redshift Lya absorbers, aris- 
ing in non-equilibrium large-scale structures that 



are highly photoionized by the metagalactic UV 
flux, as suggested by DHKW. Uncertainties from 
our current simulation limit our ability to ex- 
tract all possible quantitative information con- 
tained even in this small sample of STIS spectra, 
so improved simulations will be required to con- 
firm the results presented here. Still, these pre- 
liminary results indicate that we are close to un- 
derstanding the low redshift Lya forest at a level 
similar to what has been achieved at high red- 
shift. This would be a remarkable achievement for 
studies of the nearby intergalactic medium using 
ultraviolet spectroscopy, and analogous to what 
has transpired recently with high-z forest studies, 
promises to open up a new realm for observational 
tests of theories of cosmology, structure formation, 
and the evolution of intergalactic baryons. 
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